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ABSTRACT 

Over half of the census of massive galaxies at z~2 are dominated by quiescent stellar populations. 
The formation mechanism for these galaxies is still under debate, with models relying either on 
massive and early mergers, or cold accretion. It is therefore imperative to understand in detail the 
properties of these galaxies. We present here a detailed analysis of the star formation history (SFH) 
of FW4871, a massive galaxy at z=1.893 ± 0.002. We compare rest-frame optical and NUV slitless 
grism spectra from the Hubble Space Telescope with a large set of composite stellar populations to 
constrain the underlying star formation history. Even though the morphology features prominent 
tidal tails, indicative of a recent merger, there is no sign of on-going star formation within an aperture 
encircling one effective radius, which corresponds to a physical extent of 2.6 kpc. A model assuming 
truncation of an otherwise constant SFH gives a formation epoch zp ~10 with a truncation after 
2.7 Gyr, giving a mass-weighted age of 1.5 Gyr and a stellar mass of 0.8-3xl0 11 M Q (the intervals 
representing the output from different population synthesis models), implying star formation rates of 
30-110M Q /yr. A more complex model including a recent burst of star formation places the age of 
the youngest component at 145^7Q°Myr, with a mass contribution lower than 20%, and a maximum 
amount of dust reddening of E(B-V) <0.4 mag (95% confidence levels) . This low level of dust reddening 
is consistent with the low emission observed at 24/xm, corresponding to rest-frame 8/im, where PAH 
emission should contribute significantly if a strong formation episode were present. The colour profile 
of FW4871 does not suggest a significant radial trend in the properties of the stellar populations out 
to 3R C . We suggest that the recent merger that formed FW4871 is responsible for the quenching of 
its star formation. 

Subject headings: Galaxies: elliptical and lenticular, cD - Galaxies: individual: FW4871 - Galaxies: 
stellar content 



1. INTRODUCTION 

The formation of the most massive galaxies remains 
one of the main challenges of galaxy formation mod- 
els. Within th e observed bimodality o f the galaxy 
popu lation (e.g. iKauffmann et al.l 120031: iBaldrv et al.l 
2004), massive galaxies dominate the so-called red se- 
quence, which feature mostly old, passively evolving 
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stellar populations, suggesting an early, short-lived pe- 
riod of star formation. The implied rapid and ef- 
ficient star formation history is corrob orated by the 
presence of massive galaxies at z~l-3 (|Cimatti et alJ 
I2008t iPerez-Gonzalez et al J 12008ft and by the mild evo- 



lution with r edshift of the corresponding comoving num- 
ber density (iFontana et alJ l2006t iConselice et all [20071: 
iFerreras et all 120091: iBanerii et all 120101 ) . On the other 
hand, the strong evol ution with r e dshift found on the 
mass - size plane (e.g. [Daddi et all 120051: ITruiillo et all 
120061: lLonghetti et al.ir2007t Ivan Dokkum et al.ll2008ft re- 
mains an open question, with ongoing arguments de- 
bating whether the main size growth scenario at z< 2 
is due to either the emergence of early-type galaxies 
with cosmic time, or to intrinsic mechanisms involv- 
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observed size-mass evolution of massive galaxies has re- 
sulted in a majo r revision of our understanding of galaxy 
growth (see e.g. iNair et al.|[20TTh . 

It is believed that two main channels of galaxy 
growth are important under the hierarchical formation 
paradigm: one driven by shock heating and subsequent 
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cooling of the gas in halos (|White fc Reesl[l978l) and an- 
other one driven by the accretion of co ld gas along narrow 
filaments that does no t shock heat (iDekel &: Birnboiml 
120061 iDekel et all 120091: iKeres et all 12009( 1. At redshifts 
z>2, the presence of massive galaxies already in place 
suggests that a very efficient mechanism must operate to 
fuel and transform vast amounts of gas into stars within 
a dynamical time. 

Wide, nearby surveys along with deep surveys com- 
plemented by high resolution imaging are opening this 
important issue to scrutiny. In the past, the proper- 
ties of the stellar populations of z>l galaxies relied on 
photometric measurements or low SNR s pectra, mainly 
targe ting the rest-frame NUV (see e.g. Spinrad et al 
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1997), with the inherent degeneracies (Ferreras & Yi 



2004( 1. The advent of the Wide Field Camera 3 on board 



the Hubble Space Telescope has opened up an important 
spectral window. With its IR slitless grisms (G102 and 
GI4I) it is possible to target the Balmer break region of 
massive galaxies over a redshift range (1.5<z<3) where 
our standard ideas of galaxy growth can be put to the 
test. 

We present here a detailed analysis of the star for- 
mation history of a massive galaxy at z~2, discov- 
ered in the HST/WFC3 Ear ly Release Science by 
Ivan Dokkum fc Braiiiiueil (12010;. The slitless grism data 
from WFC3-IR has allowed us to probe in detail the star 
formation process of a massive galaxy caught in its fi- 
nal merging stages towards the red sequence. The age of 
the Universe at that redshift (around 3.4 Gyr) enables 
us to constrain in more detail the age distribution. We 
explore a large range of plausible star formation histo- 
ries, with the aim of estimating an upper bound to the 
recent and sustained star formation rate for this galaxy, 
an issue that will shed light on the main processes con- 
tri buting to the growth of massive galaxies. In contrast 
to Ivan Dokkum fc Brammcr] (|2010ft , our paper uses the 
latest flux calibration available for the WFC3 grisms, a 
very important issue when extracting star formation his- 
tories from spectral fitting. Furthermore, we explore a 
huge volume of parameter space, instead of a few trial 
SFHs, to determine not only the best fits, but robust un- 
certainties. In addition, we compare two different sets of 
population synthesis models to assess the systematics. 

In this paper we assume a standard ACDM cosmology 
with Ho=70km/s/Mpc and il m = 0.3. All magnitudes 
are given in the AB photometric system. 

2. LOOKING FOR PASSIVE GALAXIES WITH A SLITLESS 

GRISM 

Following up our previous work on the spectral analy- 
sis of massive galaxies at moderate-to-high redshift, we 
searched the WFC3/ERS slitless grism dataset, looking 
for galaxies with a prominent Balmer break. This tech- 
nique has allowed us to constrain the star formation his- 
tory of early-type ga laxies out to z^l.2 in the HUDF 
(jPascmah et al.ll2006l ) and the wider HST /ACS coverage 
of the GOODS North and South fields (jFerreras et al.l 
2009). TheIR coverage of WFC3 enables us to probe the 
very sensitive region around the Balmer break at higher 
redshifts. 

We searched through the data available from the Early 
Release Science programme (ERS) for Wide Field Cam- 
era 3 (WFC3), ID 1 1359 (PI O'Connell). A field centred 
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Fig. 1. — Slitless grism spectra of galaxy FW4871 from the Ad- 
vanced Camera for Surveys (left) and the Wide Field Camera 3-IR 
(right) . Two extractions are used in the analysis (wider extractions 
would include flux from a neighbouring object). The three panels, 
from top to bottom, show the flux, including error bars for the 0.38 
arcsec extraction, SNR and flux ratio between the extractions. The 
vertical shaded regions give the fitting interval for the analysis of 
the stellar populations. 



at RA = 03:32:17.6 and DEC = -27:42:32.4 (J2000) was 
imaged with an exposure time of 5017.6 s in each of the 
WFC3/IR F098M, F125W and F160W filters. The data 
were reduced with the standard WFC3 pipeline. Slitless 
spectroscopy of the same field was obtained through the 
WFC3/IR G102 (J band) and G141 (H band) grisms, 
for a total exposure time of 4211.7 s per grism. The 
grism spectra w ere reduced with the aXe 2.1 softwarj^l 
(Ku mmel et al.ll2010h . The data cover a field of view of 
4.65 arcmin 2 (see iStraughn et ail 120 111 for details) with 
a spectral coverage of 0.8-1.6/im and resolution between 
R=210 at A = 1/xm and R=130 at 1.4/mi. The spec- 
tral coverage allows us to probe the Balmer break in the 
redshift window l<z<3, although the shallowness of the 
ERS data restricts this range to z < 2 for the case of 
massive (M s > 10 11 Mq) galaxies, where M s is the stel- 
lar mass. 

From this dataset, we found only one galaxy with a 
high enough S/N for the application of our methodology 
to extract star format ion histories. This galaxy has bee n 
already presented in Ivan Dokkum fc Brammerl (|2010lh 
and id entified as FW487 1 from the FIREWORKS cata- 
logue dWuvts et al.ll2008l ). In this paper we look in detail 
at the constraints one can impose on its stellar popula- 
tions, employing a vast range of star formation histories. 

To further constrain the models, we also extracted the 
observed-frame optical spectrum of FW4871 from the 
2005 observations with the Advanced Camera for Surveys 
(ACS) through its WFC/G800L grism, as part of the pro- 
gram Probing Evolution And Reionization Spectroscopi- 
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Fig. 2.— Colour montage of FW4871 using WFC3/IR images in the F098M (blue), F125W (green) and F160W (red) passbands. The 
panels illustr ate the effect of the deconvolut ion. Notice the nearby distorted galaxy located NW of FW4871 is a companion (identified 
as FW4887, van Dokkum & Brammcr 20f0). The orientation of the extraction of the spectrum is shown as two white segments in the 
left-hand panel. 



TABLE 1 

Observations (AB magnitudes, errors given at the la level) 



Redshift F160W 1 R e (") n s Aperture F098M-F160W 2 F125W-F160W 2 MgUV 2 - 3 D4000 4 



1.893±0.002 19.81± 0.07 0.31± 0.04 5.4± 1.1 0.38" 1.89±0.34 0.35±0.22 1.16±0.04 1.314±0.013 

0.64" 1.88±0.26 0.40±0.16 1.15±0.05 1.297±0.012 



Total apparent magnitude corresponding to the Sersic fit. 

2 Measured within the apertures (given as diameters). 

3 As defined in lDaddi et al.1 1 120051) . 

4 As defined in lBalogh et al.1 U9991) . 

cally (PEARS, ID 10530, PI Malhotra). The total expo- 
sure time was 15140s, and the spectrum was extracted 
and calibrated with the aXe software (K iimmel et al.l 
l2009h . 

In order to avoid contamination from a nearby source 
in the grism images, we have extracted two optical+IR 
spectra for FW4871, over apertures of three and five pix- 
els on the WFC3-IR grism data, corresponding to diam- 
eters of 0.38 and 0.64 arcsec, respectively. The IR G102 
spectra have been degraded to the same spectral resolu- 
tion of those acquired with G141 (R = 130 at 1.4 /Ltm). 
The original spectral resolution of the ACS spectra is 
preserved, since these data will be fitted independently 
of the IR data. Fig. [J shows both ACS and WFC3 spec- 
tra within these two apertures. The S/N (mid panels) 
is high enough for a detailed analysis of the stellar pop- 
ulations (§4). Notice the lack of strong emission lines 
([O II], Hf3, and [O III] fall within the wavelength cov- 
erage of the WFC3 data) and the presence of prominent 
Balmer lines, typical of stellar populations that under- 
went a recent period of star formation, i.e. a k+a galaxy 
([Dressier fc Gundfl992Th 

The spectral fitting method requires knowledge of the 
effective spectral resolution. We matched a set of simple 
stellar populations with similar absorption lines (around 
an age of 1 Gyr) against the data, and smoothed the 
spectra using a Gaussian kernel with variable width. We 



also explored a range of redshifts for an accurate esti- 
m ate. We refined the reds h ift es timate from z=1.902 
in Ivan Dokkum fc"f3 rammcrl (|2010ft to z=1.893, and ob- 
tained an effective resolution FWHM of 120A, in agree- 
ment with the expected resolution of G141 (R=130 at 
1.4/zm). We note that the spectral resolution of slitlcss 
grism spectroscopy is strongly affected by object size. 
The galaxy presented in this paper is, of course, resolved, 
which implies a lower spectral resolution than the camera 
specification values, given for unresolved sources. 

3. INTRINSIC COLOUR PROFILE 

In addition to the spectroscopic da ta, we make use 
of th e NIR images from WFC3/ERS (jWindhorst et al.l 
2011). Imaging through three passbands are available: 
F098M, F125W and F160W, reaching a 50% point- 
source completeness limits of AB(F098M)=27.2 mag, 
AB(F125W)= 27.5 mag, and AB(F16 0W)=27.2 mag, at 
the 5ct level ([Windhorst et al.l 1201 ID . At the redshift 
of our target, these three passbands closely map rest- 
frame U, B and V standard filters. We ran a set of 
simple stellar popu l ations from the latest versions of the 
iBruzual fe Chariot! (|2003| ) models to determine the K- 
corrections relating these observed and rest- frame filters. 
As expected, the K-corrections are rather low, always be- 
low 0.05 mag even if the metallicity is unknown over one 
decade. 
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R/R. 

Fig. 3. — Colour profile after a Voronoi tessellation of the 
WFC3/IR images. The tessellation is performed with a target S/N 
per bin in the F160W passband above 50. The right axes give the 
K-corrected, rest-frame U-V (top) and B-V colours. 

In order to constrain the su rface brightness p rofile of 
the galaxy, we ran GALFIT ([Peng et all 120101 ) on the 
three WFC3/IR images, for a Sersic profile plus a con- 
stant sky background. The PSF is built from a set 
of 11 stars extracted from the same images. The re- 
sults are given in Tab. [TJ where the uncertainties are 
quoted at the la level and are computed from a set of 
GALFIT runs on the same images, replacing the com- 
bined image of the PSF by those from the individual 
stars. The half-light radius (0.31±0.04 arcsec) maps into 
a projected physical distance of 2.64±0.35 kpc at the 
rcdshift of FW4871, which makes this object a typical 
example of the c ompact massive ga l axies found at simi- 
lar redshift (e.g. IDaddi et al.| [200a iTruiillo et ail [20061: 
Ivan Dokkum et al. I 120081) . although we note that the 
CM05 models would locate this galaxy closer to the lo- 
cal sample on the mass-size plane (see below). The two 
extractions of the spectra we will use for the analysis 
roughly correspond to apertures with diameter 0.38", 
and 0.64", hence enclosing R e /2 and R e , respectively. 
The high value of the Sersic index makes the surface 
brightness profile compatible with an early-type mor- 
phology. 

In order to determine the intrinsic colour gradi- 
ent of FW4871, we follow the approach described in 
iFerreras et al.l (|2005l) involving a Voronoi tessellation to 
increase the S/N of the image, preserving as much as pos- 
sible the spatial information. Fig. [2] shows the original 
frames (left), and a deconvolv ed reconstru ction (right), 
using the lucy task in IRAF (jLucvl 119741 ) . The figure 
uses all three NIR passbands to generate an RGB colour 
image. This galaxy reveals itself as a major merging 
event, with nearby objects displaying tidal features as 
well, most notably the galaxy to the NW, at the same 
redshift, identified as FW4887, a galaxy with a mass 



—1/10 of FW4871, which will likely contribute to a later 
mino r merger with FW4871 ([van Dokkum fc Brammerl 
2010). The spiral-like shape of the tails might alterna- 
tively suggest a late- type galaxy morphology. However, 
as we will show below, the lack of dust and ongoing star 
formation from the analysis of the spectra makes this 
option unlikely. 

The Voronoi tessellation targeted a S/N of 50 in 
F160W, including all pixels within 3R . For instance, 
in the F098M-F160W colour profile, out of 1093 pix- 
els we get 128 binned tiles. Fig. [3] shows both colour 
profiles from the WFC3/IR data, with the correspond- 
ing rest-frame U-V and B-V colours on the right axes 
(the K-corrections were computed for a range of simple 
stellar populations at solar metallicity). We note that 
an analysis of the deconvolved images give no signifi- 
cant difference, confirming that the change with respect 
to the PSF for the three WFC3/IR passbands used, in- 
troduces no significant biases in the colors. Within the 
core (R<R e ) the colours are slightly bluer than in the 
outer regions, but still around (U-V)o ~1, reflecting ei- 
ther small levels of ongoing star formation or dust red- 
dening. The following section presents the spectroscopic 
analysis, that provides more stringent constraints on the 
underlying stellar populations. 

4. MODELLING THE STAR FORMATION HISTORY 

The high S/N of the WFC3 IR slitless grism data al- 
lows us to perform the fitting of the spectral energy dis- 
tribution corresponding to the rest-frame 3500-5500A, 
i.e. a very sensitive region to the age and metallicity dis- 
tribution. We note that a contaminant spike is present 
in the orig inal WFC3 SED at an observed A=15100A. 
At the redshift of the galaxy, this contamination appears 
very close to the metallicity-sensitive Mgl region around 
5170A (see fig. 3 of Ivan Dokkum fc Brammei1[2010l ). In 
order to be able to use this region, we corrected it by 
running the CST models (described below), fitting the 
spectrum over the wavelength range 3500-5100A. The 
best fit (extended to 5500A) is then used to perform a 
correction in the region of the spike - which extends over 
— 200A, redward of the Mgl feature. We tested our re- 
sults in the two cases, i.e. without any correction, fitting 
only blueward of 5100A, and with the correction, fitting 
out to 5500A, to find no significant difference in the age 
distributions, but an important suppression of the prob- 
ability distribution towards very low metallicities. 

The spectrum (shown in Fig. [1} shows prominent 
Balmer lines, characteristic of stellar populations of ages 
around 1 Gyr. We used the ACS data to measure the 
MgUV feature (a local peak in the heavily-blanketed 
NUV continuum at rest-frame 2625-2725A) and obtained 
a value of MgUV=1.16 ± 0.04, which is also indica tive of 
a similar age (see e.g. fig. 2 in IDaddi et al J 12005). Our 
aim in this section is to characterize the past star forma- 
tion history of FW4871. The discove ry paper of this sys- 
tem (van Dokkum & Brammer 2010) already presented a 
fir st approach to this problem, using the method outlined 
in IKriek et al.l ([20091 ) . In this paper we explore in more 
detail this issue, including additional information from 
the MgUV feature and determine robust uncertainties on 
the derived formation history by the use of a large num- 
ber of models. We ran four separate grids of models, as 
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Fig. 4. — The contours show the 69, 90 and 99% confidence levels for the spectral fitting of the WFC3 (rest-frame optical) and the ACS 
data (rest- frame NUV). In this case the models used are Simple Stellar Populations (i.e. a single age an d metallicity ) , including d ust as a 
free parameter. Th e inset shows the cumulative distribution of dust. We compare the 2007 version of the Bruzual & Chariot (2003) models 
with the models of Maraston (2005), as labelled. Two SED extractions of the WFC3 data are considered: a 3 pixel one (Left, labelled 
R<0.38", corresponding to a physical extent of 3.2 kpc) and a 5 pixel extraction (Right, labelled R<0.64", corresponding to 5.4 kpc). 



outlined below. For each choice of the star formation his- 
tory, we define a likelihood combining - as independent 
observables - the WFC3 spectral data, the measured 
colours (F098M-F1609W and F125W-F160W), and the 
MgUV index obtained from the ACS spectrum. Notice 
we have two separate runs for a WFC3 extraction of 3 
and 5 pixels, corresponding to diameters of 0.38 and 0.64 
arcsec, respectively. The ACS and photometric data are 
extracted to match the same apertures (see Tab. [T] for 
the values used) . The best fits obtained for these models 
are shown in Tab. [3] 

4.1. Simple Stellar Populations 

A first approach would involve simple stellar popula- 
tions (SSP), defined by a single age and metallicity. We 
also include dust as a free parameter described by E(B- 
V), modifying the spectrum with two reddening laws: 
the iCalzettil J2001) law - appropriate for starbursting 
syste ms - and a Milky Way reddening law (jFitzpatrickl 
1999). The grid of model parameters is defined in Tab.|H 
Fig. 2] shows the confidence levels on age and metallic- 
ity, marginalized over dust, for two population synthe- 
sis mo dels: the 2007 ver sion of the iBruzual &: Chariot! 
(120031) models (BC07; iBruzuall |2007j), and those of 
iMarastonl (|2005l hereafter CM05). Each panel in Fig. H 
corresponds to one of the two extractions, as defined 
above (other results are shown in Tab. [5J with the un- 
certainties given at the 5 and 95 percentiles). Note the 
significant difference with respect to stellar mass be- 
tween BC07 and CM05 models, a well-known feature 
caused by the different prescriptions included in these 
model s, most notably the modelling of the TP-AGB 
phase (Marast orTet al.l l2006). The insets show the cumu- 
lative probability distribution of dust, which shows good 
agreement between both models, with values below E(B- 
V)=0.02 at the 95% confidence level. Given the merger- 



like morphology of the object (Fig. [2]) it is quite remark- 
able to obtain such a low level of dust and a relatively 
"old" SSP-age (the presence of ongoing star formation 
would shift the derived age to lower values). However, it 
is a well-known fact that SSP-based fits are not robust 
indicators of age unless we are dealing with a system 
formed over a short period of time compared to stellar 
evolution timescales (e.g. a globular cluster), or for sys- 
tems with old and early-truncated star for mation histo- 
ries ( e.g. a massive early-type galaxy, see iRogers et al.l 
2010). Notice that the best fit SSP' gives a rather low 
metallicity, especially for the CM05 models. The con- 
straint on a single age to explain a composite population 
can result in such values on the metallicity. To overcome 
this, we consider more realistic composite models where 
the star formation history is parameterised by a small 
number of variables. A different reddening law does not 
change the outcome, a result that will be consistent for 
the composite models. The reason is that the reddening 
law used in this analysis pertains mainly to the optical 
region, where most reddening laws do not differ much. 
In the NUV, the only information used - the MgUV in- 
dex - is rather immune to a change in reddening law 
given that it is defined over a small spectral window, 
and it is located aw ay from the "2175A bump" (see e.g. 
IConrov et al"1l2010h . 

4.2. Exponential Models 

The lack of signatures of ongoing star formation in the 
SED of FW4871 argues against an exponentially decay- 
ing star formation rate. However, these models have 
been commonly applied to evolved stellar populations, 
and so, we present as an exercise the results of fitting 
such models (hereafter EXP models), whereby four pa- 
rameters describe a star formation history, namely the 
formation epoch (tFORi the cosmic epoch when star for- 
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Fig. 5. — As in Fig. [4] for the fit from the truncated models (CST, see text for details). The panels on the left (right) correspond to the 
BC07 (CM05) models. The top panels also show the contours for an mass-weighted average age from 0.5 Gyr to 2.5 Gyr (thick line) in 
steps of 0.5 Gyr. The composite models are based on the population synthesis models of BG07 (Jeft) and CM05 (right). 




Fig. 6. — Contours are shown at the 69, 90 and 99% confidence levels for the fit corresponding to the truncated models plus a later 
burst in the form of a simple stellar population (CST2, see text for details). The figure shows the properties of the secondary (younger) 
component, as labelled. The older component features similar constraints to the CST models shown above (i.e. Fig.[5j. In the top panels 
we show the constraints for the independent dust reddening parameters of the old (dashed lines) and the young (solid lines) components. 
The composite models are based on the population synthesis models of BC07 (ieft) and CM05 (right). 



mation starts), the timescale of the exponential (r), the 
metallicity (Z, kept fixed for all times within a model), 
and the dust content, parameterized by the colour ex- 
cess E(B-V). These so-called r-models improve over the 
SSPs as they consist of a mixture of stellar ages, giving 
a more realistic representation of the real star formation 
history, with t controlling the width of the age distri- 
bution. Tab. E] shows the range explored for these pa- 
rameters. Tab. [5] lists the parameter fits, including the 
average age and the RMS of the age distribution (using 



the likelihood as the probability distribution function of 
the parameter considered). The average mass- weighted 
stellar age according to EXP models lies between 0.5 and 
0.7 Gyr, depending on the population synthesis library 
used. We note the age of the Universe at the redshift of 
FW4871 is 3.4 Gyr. The main difference between BC07 
and CM05 models is the lower value of the age for the 
latter. Nevertheless, both population synthesis models 
agree with a rather unphysical short-lived burst of for- 
mation (t < 0.2 Gyr). In constrast to the SSP runs, 
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TABLE 2 

Parameter fits, combining the WFC3/IR SED, the MgUV index (measured from the ACS SED) 
and F098M, F125W and F160W photometry. Ages are defined as mass-weighted values. For the 

EXP MODELS. THE AGE DISPERSION At IS DEFINED AS THE RMS OF THE DISTRIBUTION, FOR CST MODELS 
IT REPRESENTS THE DURATION OF THE BURST. ERROR BARS ARE QUOTED AT THE 5-95 PERCENTILES. 



Model 
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logM s /M 
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Milkv Way Reddening (Fitzpatrick 19991 
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Fig. 7. — Best fit of the WFC3 slitless grism data for the CST2 models shown in the rest-frame, with each panel showing a different 
extraction, as labelled. The residuals (as a fraction of the RMS per wavelength bin) are shown in the bottom panel. The black (grey) lines 
correspond to the fits using BC07 (CM05) models. 
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solar and super-solar abundances are now allowed, and 
the best-fit \ 2 values are significantly reduced. These 
models would suggest that the bulk of the stellar popu- 
lations of FW4871 were formed in an intense and very 
short-lived burst started at redshift zp ~2.5-3. 

4.3. Truncated Models 

There is an important limitation with exponentially 
decaying models: the tail of the exponential creates a 
small, but significant amount of very young stars. The 
flux from those stars will affect considerably the spectral 
fit. Hence, it is necessary - for spectra without a signif- 
icant amount of star formation - to adopt an additional 
set of models, namely those where the star formation is 
truncated. For simplicity, we assume that the star forma- 
tion rate is constant within an interval between cosmic 
time tFOR and tpoR + At, and then switches abruptly to 
zero (hereafter, we refer to this set as CST models). This 
parametrisation allows us to get a more realistic estimate 
of the duration of the burst if no young stars are found 
in the spectra (as suggested, e.g. by the colour distribu- 
tion, the SSP fits of the spectra, or the absence of strong 
emission lines). Dust is also included as a free parameter 
in the same way as for the SSP and EXP models. Once 
again, Tab. [4] shows the range explored for these mod- 
els. Fig. [5] shows the confidence levels for these models, 
using the same observational constraints as for the EXP 
and the SSP models. The figure gives the constraints on 
the parameters controlling the age (top) and metallicity 
(bottom). The slanted lines in the top panels represent 
an average, mass-weighted age between 0.5 and 2.5 Gyr 
(thick line), in steps of 0.5 Gyr. Notice the degeneracy 
between age and burst duration is reduced for the anal- 
ysis with BC07 models (left panels). The CM05 models 
allow for the presence of younger stars, but the \ 2 are 
slightly worse than for BC07-based models. 

The figure suggests an average age about 1.5 Gyr, with 
a wide age spread, between 2-3 Gyr, and metallicities 
above Z /3. The results do not vary significantly with 
respect to the choice of the 0.38"(~R e /2) or 0.64"(~R o ) 
diameter apertures considered, suggesting rather homo- 
geneous stellar populations. 

4.4. Beyond Truncated Models 

One could argue that the single dust screen, adopted 
for our models so far, may not take into account a more 
complex star formation history. Complex models come 
at the price of huge volumes of parameter space. Further- 
more, the inherent degeneracies present in any analysis 
of stellar populations from unresolved spectroscopic data 
lead us to follow the approach of mapping as much as pos- 
sible the full volume of parameter space. Hence, in order 
to constrain the presence of a younger population, with 
different dust reddening to the rest of galaxy, we define 
a new set of models (CST2) consisting of the previous 
one, i.e. constant star formation truncated after some 
time lapse At, along with an extra component made up 
of a young, simple stellar population. In order to make 
the parameter search manageable, we assume a formation 
epoch for the old component of zfor = 10 and assume the 
same metallicity for both components (although allowed 
to be a free parameter). Hence, this model is described 
by six free parameters (see Tab. [4]). The new parameters 



describe the age of the young component (ty), its con- 
tribution in mass (fy) and the amount of dust in the old 
and young components (Eo and Ey). Given the small 
differences found with respec t to attenu a tion l aws, we 
will only consider a standard iFitzpatrick (|1999f ) law for 
CST2 models. 

Fig. [5] shows the constraints on the parameters mainly 
related to the young component (we do not find signif- 
icant differences for the parameters controlling the old 
component, i.e. the equivalent of Fig. [5] for the CST2 
models would be very similar). On the top panels of 
Fig. [6j the amount of dust reddening is shown both for 
the old (thin lines) and the young component (thick 
lines). Tab. [3] shows the parameter fits over the 5-95 
percentiles. Note the minimum value of the \ 2 decreases 
- although we must emphasize that all models explored 
here give acceptable fits, so that a Bayesian evidence 
criterion would not be capable of discriminating among 
these models. Nevertheless, the CST2 models show that 
a young component could be present, although with an 
age of 145^7q° Myr, and a maximum amount of dust 
reddening of E(B-V)<0.4mag (95% confidence levels). 
Its contribution to the total stellar mass budget changes 
slightly with respect to the aperture or the population 
synthesis models used, but stays roughly below 20% for 
all models, with a best-fit value in the range 5-15%. 

Fig. [7] shows the comparison between the observed 
and the best fit WFC3 spectra for the BC07 and CM05 
models, as black and grey lines, respectively. As re- 
gards emission lines, [O II] is only a lcr residual in the 
0.38 arcsec extraction, increasing to 2a in the larger aper- 
ture. H/3 is more prominent, at the 2-2. 5er level in the 
CM05 case. However, we note that FW 4871 has an X- 
ray source counterpart (|Luo et al.ll2008D and the implied 
X-ray luminosities could not be expected from a star- 
bursting system unless rates as high as l,OOOM0/yr are 
sustained, le aving us with the option of an a ctive galac- 
tic nucleus (jvan Dokkum fc Brammer1l2010h in an oth- 
erwise gas-poor merger. Alternatively, one could con- 
sider the X-ray emission originating from a hot, diffuse 
halo. This X-ray halo can be produced by the conver- 
sion of gravitational potential energy into heat. As this 
galaxy went through a major growth phase in its past, it 
is likely that this process was accompanied by a major 
growth phase via infall of material, providing gravita- 
tional potential energy. In this way, star formation could 
have continued until the h alo was heated up t o suffi - 
ciently high temperatures. iKhochfar fc Ostrikerl ()2008[ ) 
find that the fastest growing halos would also heat the 
most. In this case, the residuals found could indeed be 
related to some (low-level) ongoing star formation (al- 
though note the lack of blue colours in Fig. [3]). During 
the runs we also compute the absolute luminosity and 
the U-V colour for each star formation history. It is 
worth noting that the best fit models consistently give 
a luminosity of My = —24.89 ± 0.07 and a colour U- 
V= 1.3 ± 0.06 (AB, 90% confidence levels) regardless of 
the model or base population synthesis. 

As a final check of our (mainly spectroscopic) fits with 
respect to broadband photometry, covering a wider range 
of wavelengths, we show in Fig. [8] a comparison be- 
tween the observed photometric data and the best fit 
CST2 models obtained from spectral fitting. The obser- 



The formation of massive galaxies: FW4871 



9 



TABLE 3 

Parameter fits for the CST2 models, combining the WFC3/IR SED, the MgUV index (measured from the ACS 
SED) and F098M, F125W and F160W photometry. For ease of comparison with CST models, (Age) is defined 

HERE AS THE AVERAGE AGE OF THE OLDER COMPONENT. ERROR BARS ARE QUOTED AT THE 5-95 PERCENTILES. 



PopSyn Aperture (Age) (Gyr) At (Gyr) log (ty/Gyr) fy Ey(B-V) log Z/Z logM ij /M 



v-' 

Art 



BC07 0.38" l-55to;ig 2.37±°^ -0.85±o.3o 0.06^°^ <0.21 +0.22±°;°* 11.471°;^ 0.58 

BC07 0.64" 1-76±°;|° 1.99i°;p -0.83±°° 9 5 0.06l°;^ <0.28 +0.2lj:°;°| 11.471°;^ 0.65 

CM05 0.38" 1-86±°;|^ 1.15±^'J| -0.82t°;^ O.istjj;™ < - 10 +°- 15 ±ai6 la89 ^o:26 °- 97 

CM05 0.64" l-94±o:2? 1 - 72 -oH -0-82±°;°* 0.14±°;°;j <0.10 +0.181°;^ 10.90t°;^ 1.12 



TABLE 4 

List of parameters used in the models explored in this paper. A formation time or 

EPOCH CAN BE GIVEN EITHER BY A COSMOLOGICAL TIME (t f qr), OR BY THE CORRESPONDING 
REDSHIFT (Zpop). T(y AND Tio REPRESENT THE AGE OF THE UNIVERSE AT THE REDSHIFT OF 
FW4871 (z obs = 1.893), AND THE COSMOLOGICAL TIME AT REDSHIFT Z=10, RESPECTIVELY. 



Model 


Parameter 




Range 


Description 


No. Models 


SSP 


z FORitFOR 

log(Z/Z ) 
E(B-V) 




[10,z obs 
[-1.5, +0.3 
[0,0.5 


Formation Epoch (Single Age) 

Metallicity 

Dust reddening 


262,144 


EXP 


ZFOR;tFOR 
log(r/Gyr) 
log(Z/Z ) 
E(B-V) 




[10,z obs ] Formation Epoch 
[— 1,+0.6] Gyr Formation Timescale 
[-1.5, +0.3] Metallicity 

[0,0.5] Dust reddening 


1,048,576 


CST 


ZFORitFOR 
At 

log(Z/Z ) 
E(B-V) 


[0,(tu 


[10,z o b s ] Formation Epoch 
— *lo) ~ 3] Gyr Burst Duration 
[-1.5, +0.3] Metallicity 

[0,0.5] Dust reddening 


1,048,576 


CST2 


ZFOR, z FOR 
At 

E (B-V) 

ty 

Ey(B-V) 

fy 

log(Z/Z ) 


[0,(tu 


1C 

-tie) ~3] Gyi 
[0,0.5 
[0.1,2] Gyi 
[0,0.5 
[0,1 

[-1.5, +0.3 


Formation Epoch (Fixed) 
Burst Duration 
Dust reddening (Old) 
Young Component (SSP) 
Dust reddening (Young) 
Young Stellar Mass Fraction 
Metallicity (Both Components) 


2,985,984 



vations are given in the figure as a solid line, whereas the 
model fits are shown as open circles (BC07) or crosses 
(CM05), with error bars representing the 90% confi- 
dence levels. The observations cover a wide range, from 
the optical (HST/ACS) to the IR (Spitzer/IRAC), in- 
cluding NIR (VLT/ISAAC and HST/WFC3-IR). Notice 
the constraint from spectral fitting only use data out to 
-^obs ~ 1-6/im, whereas the photometric SED shown here 
as a solid line, comes from the total fluxes (ta ken from 
the FIREWORKS survey. IWuvts et all [2001 . In the 
case of the ACS and WFC3 data, we performed aperture 
photometry on the images, so that those fluxes corre- 
spond to the chosen 0.38 and 0.64"apertures. The IRAC 
fluxes cannot be treated this way, given the poorer spa- 
tial resolution. The fit is good for most of the photo- 
metric data points with BC07 models. In contrast, the 
CM05-based models give slightly higher NUV fluxes and 
lower 8/mi fluxes, an aspect that cannot be blamed on 
dust, since the models give a strong constraint in E(B- 
V) (see Tabs. [2] and [3]). An increase in the amount of 
dust would drastically worsen the SED fitting of the 



WFC3 data, whose superb flux calibration cannot ac- 
commodate larger amounts of reddening. The relatively 
weak flux observed in the 24/mi Spitzer/MIPS pass- 
band ($24 = 1 0-05 AtJyi compared with an 8/xm flux of 
$ 8 = 62.0 fj,Jy, IWuvts et aLl l2008D confirms that emis- 
sion from dust - and thus significant star formation - is 
rather low in this galaxy. 

5. DISCUSSION AND CONCLUSIONS 

FW4871 constitutes a typical case of a massive galaxy 
already in place by redshift z^2. This type of galaxy is 
expected to evolve into the massive early-type galaxies 
we see at low redshif t . Taking the redshift catalogue from 
the E SQ database (|Popesso et a l. 2009; iBalestra et al.l 
120101) we find 13 (54) galaxies within a redshift range 
1.8<z<2.0 at a projected separation below ±1 (2) Mpc. 
Hence, regarding environment - over the large scales 
mapped by this redshift range - this galaxy clearly sits 
in an overdense region. 

We find that a realistic model consisting of a trun- 
cated SFH gives no significant star formation in the re- 
cent past history of FW4871, whereas the duration of 
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Fig. 8. — Photometric SED corresponding to the results from the 
spectroscopic fits. The points represent the probability weighted 
fluxes according to the CST2 models, with error bars extend- 
ing over the 90% confidence interval. The line shows the pho- 
tometric data from HST/ACS, VLT/ISAAC, HST/WFC3-IR and 
Spitzer/IRAC. The solid dot with an arrow marks the observed 
frame 24/xm flux from Spitzer/MIPS (taken from IWuvts et alj 
[2001) . 

the bursting episode extends for around 2.5 Gyr. Un- 
fortunately the best fit values of the likelihood (i.e. the 
"minimum x 2 ") is similar among sets of models, prevent- 
ing us to discriminate among SSP, EXP, CST, CST2, or 
other sets of models using, for instance, a Bayesian evi- 
dence estimator, ev en with better spectral resolution (see 
iRogers et al.ll2010l for a detailed analysis in low-redshift 
galaxies) . One can only resort to prior information in the 
sense of rejecting models such as an SSP, where the for- 
mation takes place over a timescale much shorter than 
the dynamical timescales of the systems concerned, or 
EXP models, where a protracted tail of star formation 
is at odds with the lack of line emission related to ongo- 
ing star formation. Nevertheless, all models point to the 
fact that no significant star formation took place during 
the latest stages of FW4871, at least within one effective 
radius. Furthermore, the colour profile of this galaxy 
(Fig. [3]) suggests a similar outcome in the external re- 
gions, out to 3R e . 

We emphasize that within the 4.65 arcmin 2 field of 
view of the ERS/GOODS data, FW4871 is the only 
galaxy with a strong Balmer decrement observed in the 
rcdshift range sensitive to the WFC3 grism data (i.e. 
l<z<3). Given the limit in apparent magnitude of 
the observations, this constraint is compatible with the 
comoving numbe r density of massive galaxies at z^2 
(n ~!Q- 4 Mpc- 3 . iConselice etaD 12001 . Is FW4871 an 
outlier or a typical example of a forming massive galaxy? 
Larger surveys to a comparable depth are needed to an- 
swer this question. 

The analysis of FW4871 by Ivan Dokkum fc Brammerl 
( 2010) compared the WFC3-IR spectrum with a small set 



of synthetic spectra, proposing a stellar age of 0.5 Gyr 
and Ay ~1 mag, with a tentative, although very weak, 
detection of H/3 in emission in the outer regions of the 
galaxy. Their spectral fit only presented three possi- 
ble star formation histories, without any reference to 
the range of models explored. Hence, their suggested 
scenario cannot include an estimate of uncertainty in 
the properties of the u nderlying stellar populations. In 
addition, wc note that van Dokk um fe Brammerl (|2010f ) 
extracted the data from the publicly available server 
and performed the reduction of the data while the 
latest calibration files were not available. Our team 
used the best possible calibration files, obtaining a 
sli ghtly different SED compared to the one presented 
in Ivan Dokkum fc~B rammcrl (|2010D . hence the difference 
with respect to the dust content. Even though the simple 
scenarios presented in their paper are compatible with 
our analysis regarding age, their estimate of dust cannot 
be accommodated with the spectrum reduced with the 
newest calibration data. 

We extend the analysis of Ivan Dokkum fc Brammerl 
(2010) by exploring a wide range of star formation his- 
tories, and include the colour gradients to strengthen 
the hypothesis that the star formation in FW4871 has 
been truncated in the recent past of its history. The 
main conclusion from our analysis of the stellar popu- 
lations within the effective radius is that the process of 
formation started at redshift zp ~10 and lasted for about 
^2.5 Gyr. The addition of an extra, younger, component 
to the analysis allows for up to 20% in mass of a popu- 
lation with an age of ^150 Myr, with a modest amount 
of dust reddening (E(B-V)<0.4 mag), therefore with the 
bulk of stellar populations having formed at an earlier 
phase. The long At of the models along with the lack of 
star forming features in the spectrum suggests a trunca- 
tion of the star formation over a few hundreds of Myr. 
Taking into account the stellar mass of the galaxy, we 
can infer a sustained star formation rate during its burst- 
ing phase in the range 30-110 M /yr (with the quoted 
range including the difference between the stellar masses 
according to BC07 or CM05 models). One expects the 
star- forming progenitors to be sim ilar to SCUBA galax- 
ies (see e.g. iChapman et al.l [2005h . at z>2.5. FW4871 
represents the endpoint of the strong star formation ac- 
tivity of a massive galaxy, perhaps followed by a number 
of minor mergers - such as the one expected with its 
neighbour FW4887 - that will increase its compact size 
towa r ds those of local m assive galaxies (jKhochfar fe Silkl 
I2006t iNaab et al.l 120091 ). Regarding the discrepancy in 
the stellar mass given by either BC07 or CM05 models, 
we note that this disagreement is caused by the differ- 
ence in the treatment of evolved phases of stellar evolu- 
tion, with TP-AGB stars contributing significantly more 
to the luminosity in the CM05 models, thereby reducing 
the M/L especially for the ages considered. However, we 
note that this predicted con tribution still await s confir- 
mation by observations (e.g. iZibetti et al.ll2012l ). 

Given the mass of this galaxy, one should expect a 
major contribution t o its formation via cold a ccretion at 
high redshift (see e.g. Dekel & Birnboim 2006), a process 
that would boost the star formation rate both from the 
enhanced gas accretion rate and a possible increase of th e 
star formation efficiency (see e.g. lKhochfar fc Silk 2009). 
The derived star formation history of FW4871 provides 
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evidence - via the superb spectroscopic data from HST 
- towards a highly efficient channel of galaxy formation 
that started around zp ~10, and was subsequently trun- 
cated around z>2, followed by a brief quiescent phase. 
What caused this truncation? Two main scenarios can 
be brought forward: i) The onset of a galactic super- 
wind from stellar feedback: a sustained star formation 
rate of about ~100M Q /yr for a couple of Gyr will pro- 
vide enough energy; ii) The previous merging process 
that resulted in the formation of FW4871 could be the 
ultimate caus e of the quenching, p ossibly by switching 
of the AGN (jHambrick et all 120111 ). The morphologi- 
cal appeara nce of FW4871, and the X-ray detection of 
this source (|Luo et al.ll2008l ) points to the latter. In this 
case, the recent merging process to form FW4871 could 
be considered as the cause of truncation. One could en- 
vision FW4871 as a high redshift equivalent of a post- 
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starburst, or k+a galaxy. The analysis of the environ- 
ment in a sa mple of k+a galaxies from the Sloan Digital 
Sky Survey (|Gotol l2005h reveals that the truncation of 
the star-bursting phase is caused by the interaction with 
a nearby companion, a similar mechanism to the one con- 
templated here. 

This process would make major mergers an important 
factor to explain the "road to the red sequence" , namely 
the process of quenching of star formation in massive 
galaxies, an issue that may require proper calibration of 
semi-analytic models of galaxy formation. 
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